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ABSTRACT 

We compute the one-loop effects on the neutrino propagation through matter induced by 
virtual supersymmetric particles. We show that, in the minimal version of the supersymmet- 
ric standard model, a non- degeneracy between sleptons of the second and third generations 
can have sizeable effects on the v^-Vr oscillations in matter. In particular, we discuss how 
this could affect the detection of the energetic neutrino fluxes arising from annihilation of 
supersymmetric dark matter in the center of the sun. 



The fact that ly-boson exchange with electrons in ordinary matter affects the Ve index of 
refraction but not those of the other neutrino flavours is the basis of the Mikheyev-Smirnov- 
Wolfenstein (MSW) effect which provides the nicest explanation to the solar neutrino 
problem. 

In the Standard Model (SM), at small momentum transfer muon and tau neutrinos 
interact indistinguishably with ordinary matter at the tree level. However, a difference 
among the and Vr indices of refraction in matter appears at one-loop [3, although it is 
suppressed by O 'W^ with respect to the size of the charged current effects affecting 

the Ve propagation, having then probably no observable implications. 

In this paper we compute the one-loop contributions to the neutrino refraction indices 
in the best motivated extension of the SM, i.e. the minimal supersymmetric version of it 
(MSSM). We show that there are potentially much larger radiative effects than in the SM 
itself and we discuss the possible physical relevance they could have. 

The interactions of neutrinos with matter are described by the matrix element 



For neutrinos propagating through an unpolarized medium at rest, the temporal compo- 
nent of the fermionic vector current leads to a non-vanishing neutrino forward-scattering 
amplitude and, hence, to a neutrino refraction index riy given by [BI 



f=u,d,e 



Nf is the number density of fermion / and, at the tree-level. 



T3{fL)~2Qfsl, + 6, 
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with 



sin^Ow, T-ii^fi) the third component of isospin of fi and Qj its charge. 



The indices of refraction affect the neutrino flavour evolution during propagation, which 
is described by (for reviews see [||) 
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where V is the unitary matrix relating the neutrino flavour (z/^, with a 
mass (z/j, with i = 1, 2, 3) eigenstates, i.e. z/^ = Vaii^i. Also Am?- = m 
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Auajj = n^^ —rii,^. (We note that, although the indices of refraction will be computed in the 
MSSM, some departure from it should be responsible for the neutrino masses and mixings 
themselves). 
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Since matter effects in oscillations involving Pf, will be largely dominated by the charged 
current (CC) term 5if in eq. (3), the radiative corrections to the tree-level result Arig/x = 
— \piGF^elVv turn out to be negligible. We will then concentrate on the computation of 
/^rir^i which does not involve the CC piece and vanishes at the tree- level. 

It is useful to parametrise the effects of the radiative corrections, for / 7^ ^, as 

Clf = p^^fUfL)~2QjX'^^fsl. (5) 

The splitting of the radiative effects between p and A is somewhat arbitrary, and it is con- 
venient to include in p the /-dependent box diagram contributions. This has the advantage 
that, for a neutral medium, the corrections in A (which include also contributions from the 
neutrino charge radius 0) turn out to be /-independent (see ref. 0) and do not contribute 
to Arirn, due to the fact that 

E NjQf = 0. (6) 



f=u,d,e 



One has then 



p.Aur^ = -V2Gf NfTs{fL)Apf, (7) 
/ 

with Apf = p^-f - p"''^. 

In the SM, Ap^ gets contributions from the one-loop corrections to the uuZ vertex and 
VT-boson box diagrams 0, leading to 
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»7r 

with X = m^/M^ and aw = ct/s^. These corrections are small due to the one-loop factor 
aw/ 4:71, but are also quite suppressed by the smallness of the factor m'^/Myy ~ 4 x 10~^. 

The computation of the supersymmetric contribution to Ap requires the evaluation of 
the Feynman diagrams depicted in fig. 1, leading to the results that are summarized in the 
Appendix. The important point is that, besides the one-loop factor aw/4:TT, what fixes now 
the size of Ap is mainly the splitting among the sleptons of the second and third generations. 

At this point, it is useful to recall that a usual simplifying assumption made in phe- 
nomenological applications is to consider all sfermions to be exactly degenerate at the GUT 
scale, and obtain their low energy splittings from the renormalization group evolution of 
the soft parameters and from terms arising after the electroweak symmetry breaking. In 
this way, although squarks get significantly splitted from sleptons, the splittings among the 
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masses of different slepton generations are only due to the small r- Yukawa coupling. This 
usually implies that m? — m? is O(m^), and hence the radiative effects on the z/^,,- indices of 
refraction are in this case not larger than the SM ones. An exception to this, still assuming 
a universal soft scalar mass m, is when there is a large tl-th mixing. This happens for large 
values of the Higgs mixing parameter and large tg/3, or for large values of the parameter A 
of the trilinear soft terms, in which case the splitting can be 0[mT-(Am + /itg/5)] (see ref. 0). 

From a more general perspective, the universality assumptions (which give the easiest 
way to get rid of FCNC phenomena) are not really a necessity, and actually non-universal 
soft terms usually arise in string theories and can also be generated in GUTs Universal 
sfermion masses may not even be desirable in some respects, and it has been argued that 
non- universalities may prove useful in reconciling different phenomenological constraints in 
super symmetric GUTs p|, §]. Also, it has recently been suggested that sfermion masses may 
dynamically align along the directions, in flavour space, of the fermion masses, suppressing 
FCNC but allowing large mass splittings |10[| . 

If one considers the general case in which a sizeable splitting is allowed among ^ and r 
sleptonsQ, the SUSY contribution to An^^ could then be larger than the SM one. 

We will present an illustrative situation in which the effect here described is important and 
then comment on how the results are modified when one changes the starting assumptions. 
For simplicity we assume no Jl^Jr mixings, neglect intergenerational mixings of sleptons as 
well as splittings due to D-terms among charged and neutral sleptons or among ii and Er, 
which are anyhow inessential to the conclusions reached. We take first generation sleptons 
degenerate with the second generation ones, and only allow the third generation sleptons 
to have a different mass. We take a light but experimentally allowed value for the second 
generation slepton masses, = mc,^ = Max[60 GeV, + GeV], where is the lightest 
neutralino mass (which we assume to be the lightest supersymmetric particle (LSP)), and 
take second and third generation sleptons to be splitted by an amount — = 60 GeV. 
We also assume that squarks are much heavier than sleptons (as usually results from the 
effects of gluino masses in the renormalization group evolution of scalar masses). This last 
implies that box diagrams only contribute sizeably to interactions with electrons. 

In fig. 2 we plot the ratio of the SUSY and SM values of An^^ for an isoscalar medium 

{Yn = Nn/Np = 1), as a function of the supersymmetric parameter space that determines the 

chargino and neutralino masses and couplings (SU(2) gaugino mas^ M and Higgs mixing 

parameter /x). We present results for values of the ratio of Higgs VEVs tg/5 = V2/V1 = 2 

^This splitting is not directly related to very suppressed rare processes such as /i ^ 67, ^ — ^ 3e, t ^ ^7, 
etc. but could give rise to small universality violations p^ . 

^for definiteness we assumed common gaugino masses at the GUT scale to obtain the neutralino spectra 
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(fig. 2. a) and tg/5 = 40 (fig. 2.6), showing that the dependence on it is only mild. The dark 
regions for small values of fi and M are excluded by the LEP constraint > 45 GeV, 
that is the main bound from accelerators. 

It is apparent from fig. 2 that, for the slepton mass splittings considered, Ari^^ may 
be an order of magnitude larger than in the SM. The SUSY contribution turns out to be 
dominated by the chargino boxes and penguins involving El exchange. The neutralino boxes 
are generally small, while neutralino penguins, not shown in fig. 1, give no contribution 
(similarly to what happens for instance in 6 ^ si^£~ decays [1^). Thus, the relevant 
splitting is the one among charged sleptons rather than among sneutrinos. The effect is 
especially large in a region of parameter space where the chargino masses are below 80- 
100 GeV, i.e. testable at LEPII, and becomes less important for large values of \fi\ and M, 
i.e. for heavier charginos. Slepton splittings smaller than the one adopted would lead to 
proportionally smaller effects, while larger splittings can increase the effect by up to a factor 
of two. The neglected box diagrams involving squark exchange may also increase the SUSY 
contribution. Also note that, since penguin contributions to Ap are /-independent, they 
lead (see eq. (7)) to a contribution to An,-^ proportional to Ne{2 — F„) (in the sun, F„ varies 
from ~ 0.16 in the surface to 0.5 in the center). 

Clearly the sign of An^^, and hence whether resonant matter effects take place among 
neutrinos or antineutrinosQ, depends, for significant slepton splittings, on whether f are 
heavier or lighter than fi {nif smaller than leads to a resonance crossing among neutrinos 
if Am|2 > 0). 

Let us also note that the penguin and box diagrams involving exchange are propor- 
tional to the square of the lepton Yukawa coupling, so that their effect on n^^ is much smaller 
than that on n,y^ . This fact has the interesting effect of making their contribution to An,-^ to 
depend just on rrif, rather than on a slepton mass splitting. However, even for tg/3 = 40 and 
= 60 GeV, they give a contribution not larger than the SM one. Similar conclusions hold 
for the penguins involving exchange (keeping in mind that in the MSSM mH+ > M^). 

We finally mention that other extensions of the SM may also sizeably affect Aurfi- In 
particular, supersymmetric i?-parity violating interactions can modify the neutrino indices 
of refraction already at the tree-level |T^, although those models would be less interesting 
as regards the application discussed below. Another simple example would be the presence 
of a new Z' gauge boson with non-universal couplings to leptons [[15|] . 

We turn now to consider the possible physical relevance of these radiative effects for u^- 
Vj. matter oscillations. When discussing applications, we will neglect the z/^^^ mixing with 
•^for antincutrinos, the sign of the matrix element is reversed, so that rip^ — rtp^ = — {n^^ — n^^) 
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Ve to be left with just a two flavour situation. The generalization to three flavour neutrino 
mixing should pose no problems. 

A first difficulty to observe any conversion among and Vj. is that for low energies, < 
0.1 GeV, these neutrinos are only detected by means of their neutral current interactions and 
are hence essentially indistinguishable. Furthermore, only e-type neutrinos are produced in 
the sun (except for possible ordinary MSW conversions inside the sun) and equal amounts 
of z/^ and Ur are produced in supernovae, so that oscillations among them do not give 
actually any overall result. These problems are not present in long-baseline z/^ oscillation 
experiments on earth (either with u^^ from accelerators or using atmospheric neutrinos), 
but it is easy to convince oneself that the resonance oscillation length in terrestrial matter 
(inversely proportional to Ara,-^) is typically much larger than the earth diameter, and hence 
oscillation effects are negligible. 

The situation that we want to describe, in which the matter effects here analysed are 
relevant, is actually directly related to the supersymmetric framework under consideration. 
It is well known that a nice feature of supersymmetry, once i?-parity conservation is adopted 
to avoid B and L violation, is that the LSP, usually a neutralino, is stable and naturally 
becomes a good dark matter (DM) candidate. There are two main strategies that are being 



pursued at present to experimentally search for SUSY DM |]T6|. The first is the direct 



search of the energy deposited by halo neutralinos interacting with target nuclei, and the 
second one is the search of energetic neutrinos produced in the annihilation of DM trapped 
in the interior of the sun or the earth [00. In particular, upward-going muons produced in 
the rock (or ice) just below underground detectors by energetic z/^ and (with E^, > few 
GeV) may allow to probe significant regions of the supersymmetric parameter space in new 



installations such as Superkamiokande or Amanda. As was shown in ref. |T8|, the usual 
MSW effect between and u^^r can affect the detection rate predictions. 

To show the possible effects of energetic v^-Vt niatter enhanced oscillations in the solar 
interior, we plot in fig. 3 the z/^ neutrino survival probability in the sin^2^^T- vs. Am|2 plane, 
assuming that e = AriT-^/An^g = 10~^ (in the SM, e ~ —5 x 10^^). The contours correspond 
to P{Vfj, —>■ Vfj) = 0.8 (continuous lines) and 0.45 (dashed lines) for two neutrino energies^]. 
El, = 10 and 40 GeV. It is clear that, for significant ranges of Am^ and sin^ 29, the oscillations 
of high energy neutrinos are sizeably affected by matter effects. For decreasing values of |e|, 
the MSW type resonant effects take place for smaller Am^ values. The adiabatic condition 
in the resonance transition becomes harder to achieve, making the regions of significant 
"'halo X annihilations may also provide some signals. 

^for energies above 100 GeV, i' absorption in the sun starts to be relevant, making the oscillation formalism 
to be no longer valid. Our computation of the box diagrams, which neglected external momenta, actually 
contained terms proportional to pu ■ Pf, which would modify the results only for Ei, > few TeV. 



6 



transition to shrink towards large mixings, and the effect eventually becomes very small for 
|e| < 10-^ 

These v^-v^ oscillations may have important implications for the detection of neutralino 
annihilation signals. This is because the and fluxes from neutralino annihilation are 
generally quite different, so that oscillations among them modify the expected signal 
at underground detectors. The difference among and fluxes has its origin in the 
fact that the non-relativistic neutralino annihilation cross section into fermion pairs // is 
proportional to m^, either due to a p- wave suppression (as for annihilations mediated by Z or 
sfermion exchange) or due to a Yukawa suppression (as in the case of Higgs boson mediated 
annihilations). Hence, neutralinos do not directly annihilate into neutrino pairs. Different 
neutrino fluxes result then from 6, c and r decays |]19| (for < M^y, since otherwise other 
channels involving gauge bosons in the final state are also allowed and can produce prompt 
secondary neutrinos of different flavours in similar amounts). Furthermore, light mesons and 
muons produced in the x annihilation are stopped by the solar medium before they decay, 
yielding no secondary fluxes of energetic neutrinos. 

Rather than scanning all the supersymmetric parameter space, we will consider as an 
illustrative example the simple but still quite general case in which the lightest neutralino 
is mainly gaugino, i.e. > M, with < M-\jy (in this region the effect here discussed is 
potentially large and also the neutralino cosmological relic density is usually significant). If 
squarks are heavier than sleptons, as we are assuming, and sleptons are not too heavy, the 
main non-relativistic neutralino annihilation channel is by t-channel f exchange, producing 
a rf pair. The and z/^ fluxes from the subsequent r decays will then clearly be quite 
different. (In the general case of an arbitrary neutralino composition and squark masses, 
the fluxes are still different but one needs to include the extra annihilation channels and the 
model dependent branching ratios entering in the v^.^ yields). 

In fig. 4 we show the and z/^ differential neutrino yields (with thin dashed and solid lines 
respectively) produced by the annihilation of neutralinos into a r pair (the main annihilation 
channel in our example). What is actually plotted is z^dN/dz (where z = E^/m^), which 
is the relevant quantity for underground detection^ because both the CC cross section 
and the muon range in the rock (or ice) are proportional to the neutrino energy. 

With thick lines we show how the fluxes get modified after traversing the solar interior, 
assuming = 50 GeV and taking (e, Am^, sin^ 29) to be (10~^, 6 x 10""^ eV^, 0.1) in fig. 4. a, 
(10-^ 3 X 10"^ eV^, 0.1) in fig. 4.6 and (10-^ 10"^ eV^, 0.6) in fig. 4.c. The neutrino masses 
and mixings assumed in the case of fig. 4.c lie in the region, also shown in fig. 3, of interest 
to explain the atmospheric neutrino anomaly. One should note, however, that for large 
^the muon flux arising from 1/^ interactions is quite suppressed 
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mixings the effects of vacuum oscillations are significant and would affect both neutrinos 
and antineutrinos, while matter induced oscillations affect either one or the other. 

From the results, one can see that the detection rates may be sizeably modified by the 
matter enhanced oscillations here described, providing an interesting physical manifestation 
of the radiative supersymmetric effects studied. The uncertainties involved in the theoret- 
ical predictions of the DM annihilation signal (unknown neutralino mass and composition, 
uncertainties in the local halo density and DM velocity distribution, etc.) will however com- 
plicate the interpretation of any positive detection, so that these effects should actually be 
considered as providing a further spread in the theoretical predictions until these parameters 
become more constrained by accelerators and direct DM searches. The shape of the neutrino 
spectrum gives probably a clearer signature of the matter oscillations, and may provide a 
useful handle to identify them. 

I would like to thank Stefano Bertolini for very useful discussions. 
Appendix 

We summarize here the supersymmetric contributions to An^-^ arising from the diagrams 
in fig. 1. 

A few simplifying assumptions will be adopted. Motivated by the smallness of FCNC 
phenomena, we will assume that sfermion masses align, in flavour space, in the directions 
of the fermion masses. This makes the neutralino-fermion-sfermion vertices diagonal in 
generation space. We will also ignore mixings, so that the appropriate mass eigenstates 

are fiL and fm. These assumptions simplify the calculations but are not essential to the 
conclusions reached. 

Following the notation of ref. [0, we denote by Z'^j the 4x4 matrix diagonalizing the 
neutralino states in the basis (7, Z, Hi, H2), and U and V are the 2x2 matrices required for 
the diagonalization of the chargino mass matrix. In the radiative corrections involving the 
neutralinos only the gaugino components will contribute sizeably since the higgsino cou- 
plings are very small for interactions with ordinary matter. The Feynman rule for the fx^fi 
vertex can be then parametrised as —ig\/2G-'fj^PR while that involving by ig\^G-'f*j^PL. 
The couplings to the gaugino components are given by 
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with c{(^) 



Regarding the chargino couphngs, we will ignore Cabibbo-type intergenerational mixings 
(for leptonsQ and, in the box contributions also for quarks), since the flavour-conserving 
processes under study are not GIM suppressed. We then only include the chargino mixing 
matrices in the vertices (e.g. the v^XjC-i vertex is —igUjiPR, the i^eXj'^R vertex is igHjPfi, 
with H^j = meUj2/V2Mwcp, see ref. M). 

The penguin and box contributions to Ap-^ can be written as 



box 



The (/-independent) contributions to the penguins involving charginos and II exchange, 
arise from the self-energy diagrams in fig. l.a, App(E), and from the diagrams where the 
Z couples to the slepton-line App{i) (fig. 1.6), or to the chargino-line, App(x+) (fig. l.c). 
Direct computation leads to the results: 



Ap^(S) = ^E|t^: 
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where X^h = {iTia/mh)'^ and fi is an arbitrary mass scale. The couplings O'^j = —VnV*-^ — 
\Vi2V*2 + dija^Y and O'f^ = —UiiU*-^^ — ^Ui2U*2 + SijS^y describe the xtXjZ vertex, that reads 



Cw 



Y[0%Pl + O'f^Pn]. The functions Fq and Go are 



Goix,y) = 



xlnx 



{x - y){x 
x^lnx 



\x 



y){x 



+ (x ^ y), 



y) 



3 
2' 



The penguins involving exchange, that although proportional to {mr/MwY '^ay in 
principle be enhanced for large values of tg/5, give a total contribution (neglecting the jiR 
^this may not be a good approximation for neutrino mixings close to maximal 
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exchange that is proportional to m^/M^) 

2 r 



3=1 



3 



Similarly, the charged Higgs boson contribution in fig. l.a-c is 



1 ^ Iny 



1-y {l-yY 



where y = {mr/Mff+y. 

The box diagrams (fig. l.d) involving charginos and £l exchange give: 



aw M? 2 



4% ml . , -, 

i^e J,k = l 



dr. Ji'!— 1 



dL J, 



The boxes with in exchange are just obtained by replacing in the previous expressions 
£l — ^ UjiUki — > HJHl* and omitting the flu exchange contribution. 

Finally, box diagrams involving neutralinos (fig. l.e, /) interacting with f = e, u, d give: 



ApL(x°) - — ^ E GtiGli I ^ <^/i'^/iv'^x°/L^xg/L^'(^^./L'^x°/L'^x^/J 



M2 



where 

F'(x,y,z) = 

G'{x,y,z) 



Qj Q*k 

2 ^ ^^i>rfL ' ^X?/i. ' ^Xg/i) 



x\n.x 



+ {j ^ k, L ^ R)\ - {ur Ui^), 



y\ny 



z\nz 



{x - y){x - z){x - 1) {y - x){y - z){y - 1) {z - x){z - y){z - 1)' 
x'^lnx y^lny zHnz 



+ 



+ 



(x - y)(x - z){x - 1) (y - x)(y - z){y - 1) (z - x)(z - y)(z - 1) 
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Figure Captions 

Fig. 1: Feynman diagrams describing the supersymmetric contribution to z/^-/ forward scat- 
tering {i = /i, t). The blob in fig. l.a represents corrections to both u external legs. 

Fig. 2: Ratio of the supersymmetric and SM values of Aur^ for an isoscalar medium. We 
take the sleptons of the first two generations to have a common mass of Max [60 GeV, 
rriy. + 20 GeV] (see text), and assume the third slepton generation to be heavier by 60 GeV. 
Fig. 2. a is for tg/5 = 2 while fig. 2.6 is for tgjS = 40. 

Fig. 3: Contours of survival probability P{i^^ u^) = 0.8 (solid lines) and 0.45 (dashed 
lines) for neutrino energies E^, = 10 and 40 GeV, taking e = Arirn/ Anfj_e = 10~^. 

Fig. 4: Differential (solid lines) and (dashed lines) yields (x^;^) vs. z = Eiy/m^, 
for X annihilations into r pairs. Thin lines describe the original spectra (ref. |[T^) while 
thick lines include the matter effects, assuming = 50 GeV. The figures correspond to 
(e, Am2,sin2 2^) equal to (lO'^, 6 x 10"^ eV^, 0.1) in fig. 4.a, (10-^ 3 x 10"^ eV^, 0.1) in 
fig. 4.6 and (10-^ 10"^ eV^, 0.6) in fig. 4.c. 
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